Optical conductivity in doped manganites with planar x 2 -y 2 orbital order 
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We investigate a planar model for the ferromagnetic (FM) phase of manganites, which develops 
orbital order of e g electrons with x 2 -y 2 -symmetry at low temperature. The dynamic structure factor 
of orbital excitations and the optical conductivity u(u>) are studied with help of a finite-temperature 
diagonalization method. Our calculations provide a theoretical prediction for cr(a>) for the 2D FM 
state and are of possible relevance for the recently found A-type phase of manganites at high doping 
which consists of FM layers coupled antiferromagnetically. In the x 2 -y 2 ordered regime <x(w) shows 
both a Drude peak and a gapped incoherent absorption due to a gap in the orbital excitations. 
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Recently Kawano et al. M and Akimoto et al. || ob- 
served a new ground state in heavily doped 3D mangan- 
ites Pr 1/2 Sr 1/2 Mn0 3 and (Lai_ z Nd 2 )i_ 3; Sr a; Mn03 (e.g. 
x = 0.54 and z = 0.0 — 1.0), which surprisingly show 
the same antiferromagnetic A-structure known from the 
undoped LaMnC>3 where ferromagnetic (FM) planes are 
coupled antiferromagnetically. It was proposed that the 
occupied orbitals in this phase have x 2 -y 2 symmetry [0 in 
contrast to the alternating orbital order in LaMnC>3 ||, 
and as a consequence the conduction should be quasi- 
2D and the magnetic coupling between the layers small. 
This broken symmetry state has to be distinguished from 
the orbital structure of the 3D FM state realized at lower 
doping, where orbital correlations are expected to be liq- 
uid like • The latter FM state shows a highly anoma- 
lous form of the optical conductivity |^-|| with a large 
incoherent absorption extending up to uo ~ leV. It has 
been argued that the anomalous incoherent absorption is 
related to the orbital degree of freedom @J||| and char- 
acteristic for an orbital liquid However there are 
also alternative scenarios which attribute the incoherent 
feature to Jahn- Teller polarons Comparable opti- 
cal studies in the high doping regime where the orbital 
ordered A-phase is the stable ground state have not yet 
been reported. The purpose of this work is to present a 
theoretical analysis of the possible outcome of such exper- 
iments. The frequency dependence of o~(u>) is interesting 
since it is expected to reveal the orbital excitations in 
the doped x 2 -y 2 ordered phase, and may allow to distin- 
guish between the orbital excitation and the Jahn- Teller 
polaron scenario. 

The evolution of x 2 -y 2 orbital order upon doping was 
also found in the study of planar models derived from 
the degenerate Kondo lattice model in the fully polarized 
FM phase |L0]]. Therefore we believe that this model is 
an appropriate starting point to investigate these ques- 
tions. The i2g-spins of the Mn 3+ and Mn 4+ ions as well 
the e ff -electron spins of Mn 3+ are aligned globally in the 
FM-phase and can be integrated out. What remains is 
the orbital degree of freedom since there are two nearly 



degenerate e g orbitals per site. The large intra-atomic 
repulsive interactions prevent double occupancy of e g or- 
bitals at the same site, which leads to a strong correlation 
problem and to nontrivial charge and orbital dynamics. 

For the saturated FM state the restriction to configu- 
rations without double occupancy leads to the orbital t- J 
model |l]@ 
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with d\ a = <ii a (l — n-ia)- Here we use a and b (a and 0) as 
orbital-pseudospin indices; while a denotes the orthog- 
onal e g orbital with respect to orbital a. A convenient 
basis is | f) = d x i_ y i and | |) = d 3z 2_ r 2. The transfer 
matrix elements are then given by 
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which allows for inter-orbital hopping in the icy-plane. 
The =F-sign distinguishes hopping along x and y direc- 
tion. Here £/4 is the matrix element between d 3z 2_ r 2 
orbitals in the xy-plane. For one electron per site this 
model describes a Mott insulator. The orbital interac- 
tion Hi n t follows as a consequence of the elimation of 
doubly occupied sites with energy of an effective U jl2) , 
where we assume that because of the large Hund-splitting 
the high-spin intermediate states are most important. 
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Here u, u' = (±a, 0) or (0, ±a) are lattice unit vectors. 
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Hfa t consists of 



The orbital interaction Hi nt - 
two types of contributions: (i) 2-site terms (u = u'), i.e. 
similar to the Heisenberg interaction in the standard t- 
J model, yet more complex because of the nonvanishing 
off-diagonal t ab , and (ii) 3-site hopping terms (u ^ u') 
between second nearest neighbors. In the half-filled case 
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only the 2-site interaction is operative and may induce 
some orbital order. In the presence of hole doping, i.e. 
for less than one e g -electron per site, both the kinetic 
energy and the 3-site contributions in H int lead to prop- 
agation of the holes and to a frustration of the orbital 

(2) 

order prefered by the orbital interaction H^ n ' t . When ne- 
glecting the 3-site terms our Hamiltonian is equivalent to 
the orbital model studied by Ishihara et al. jll] . Since the 
spectral shape of a(uj) is only weakly affected by these 
terms jlO| , we do not consider them in this study. For the 
calculation of correlation functions (CF) we use a gener- 
alization of the exact diagonalization technique for finite 
temperature developed by Jaklic and Prelovsek UM- 
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correlation effect. 
responding to t ^ 



Results are shown for U/t = 16, cor- 
0.3 eV and U ~ 4.5 eV [pi. 
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FIG. 1. Evolution of orbital correlations (T?T? +R ) as 
function of distance R for various temperatures T. The ratio 
of hole number and system size is Nh/N=10/16 and U = 16i. 
The inset gives a schematic view of the orbital correlations 
for T = 0.05i (strength proportional to radius). 

Figure 1 monitors the evolution of x 2 -y 2 orbital cor- 
relations as function of distance and temperature for a 
16-site cluster Q. Here we use the pseudospin operator 

de- 



notation 2? = ^(riif — riii ), = ^rf^U and T; 
There is a clear manifestation of long-range orbital cor- 
relations below T = 0.5t which saturates below T — Q.2t. 
Although {T'T'.-g) ~ 0-0^ appears small, it is close 
the maximal value ~ 0.03 for this CF given the elec- 
tron concentration 6/16. This CF has a minimum for 
nearest-neighbors (see inset of Fig.l), which is partially 
due to the evolution of an exchange-correlation hole in 
the charge CF (nini+R). In addition (nin- 1+ ji) shows fur- 
ther neighbor density modulations, which disappear with 
the orbital order. Usually such charge correlations near 
quarter-filling are attributed to strong nearest-neighbor 
(n.n.) repulsion. Since such interactions are not con- 
tained in (1), charge correlations appear here as pure 
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FIG. 2. Transverse orbital excitation spectra (T^TZq) 
(a) in the orbital ordered low-temperature regime for various 
temperatures and (b) in the high-temperature regime with 
gapless orbital excitations for same parameters as in Fig. 1. 
Spectra are broadened using r = 0.0 If. In (a) a comparison 
with results from the bandpicture is included (dash-dotted 
lines) . 

The orbital excitations are crucial for the understand- 
ing of the optical conductivity of manganites in the FM 
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phase. The form of the orbital excitation spectra changes 
significantly with doping. In the undoped case with alter- 
nating orbital order due to the interaction (3) the orbital 
excitation spectra (T^TZq) are essentially described by 
two collective modes pl|, while at low doping with the 
emergence of x 2 -y 2 orbital order this changes into a single 
collective orbital mode which is gapped. In the high dop- 
ing case, x ~ 0.6, studied here (Fig. 2a) the single mode 
features have almost disappeared and the small cluster 
results show the features of a continuum with spectral 
weight distributed up to energies of about ~ At. We 
note, however, that the dominant peaks in these spectra 
correspond to the collective excitations observed for sys- 
tems at low doping. 

In the low and moderate doping regime a slave fermion 
approach provides a satisfactory description of the or- 
bital excitations. In such a treatment holes are described 
by spinless fermions h\ and the orbital pseudospins by 
Schwinger boson operators b- la leading to the kinetic en- 
ergy Hamiltonian H t = - X)(y> a 6*if ^M<A'i> + H - c - 
Subsequent mean-field treatment, i.e. assuming con- 
densation of bit, plus Gaussian fluctuations leads to 
H ~ X^q^q^q w ith uiq = 3xt(l — 7 q /3), where 7 q = 
(cosq x + cosq y )/2. This result tells us that the energy 
scale of orbital excitations is proportional to the doping 
concentration x, which is consistent with the diagonal- 
ization data. The scaling of w q with x is due to the fact 
that x 2 -y 2 orbital order is induced by the kinetic energy of 
the holes which overcomes the superexchange interaction 
(3). Further boson- fermion coupling terms not included 
in this consideration lead to deviations from single mode 
behavior and will be discussed elsewhere. 

In this work we provide a comparison with an ana- 
lytical approach based on the band picture which was 
proposed by Shiba and coworkers ||. This Fermi- liquid 
approach assumes that U/t is not too large in sytems like 
Lai_ 2; Sr 2 ;Mn03. Moreover the approach is expected to 
be useful in the dilute (e g electron density) limit. In this 
case the constraint in (1) may be neglected and correla- 
tion effects can be included in a perturbative way. Thus 
the kinetic energy can be diagonalized, leading to two 
tight-binding (TB) bands. For the planar model the en- 
ergies are E±/t = 2jk ± + 3?y£ — /j,, where \x is the 
chemical potential and rj^ = (cos k x — cos k y )/2 . Results 
for the transverse orbital excitation spectrum (T+TZq) 
are shown in Fig. 2a for T = and compared with 
the diagonalization data. The TB-spectra for the planar 
model, i.e. neglecting correlation effects, show a broad 
continuum with several van Hove singularities. The pro- 
nounced structures in the diagonalization data around 
lj ~ 3xt, reminescent of the collective orbital waves at 
lower doping, are not reproduced. 

At temperatures T > 0.2 thermal fluctuations destroy 
the orbital order (Fig. 1). In this regime the orbital ex- 
citation gap changes into a pseudo gap with increasing 
low energy spectral weight (Fig. 2b). The appearance 



of low energy excitations in the orbital liquid phase has 
a dramatic effect on the optical conductivity as we shall 
see in the following. 

The frequency dependent conductivity consists of two 
parts [^6|, cto(w) = 2-Ke 2 D c 8(uj) + a(w). The (5-function 
contribution is proportional to the charge stiffness D Cl 
which vanishes in insulators. The finite frequency absorp- 
tion (or regular part) <j(u>) is determined by the current- 
current correlation function: 

1 - e-"/ T 
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-Re I dte^(j x (t) Jx ). (4) 



where the x-componcnt of the current operator is given 
by = -ieXj+uafc^+uj u * dt +ua d jb . 



3 
o 

o 

a. 



T/t=0.05 
T/t=0.1 

T/t=0.2 i 0.10 




0.0 1.0 2.0 3.0 4.0 5.0 



3 °- 3 

o 

a. 

0.2 



(b) 



T/t= 


0.05 


I T/t= 


0.2 


I T/t= 


0.5 


|l T/t= 


1.0 















1.0 2.0 3.0 4.0 5.0 

(0/t 




T/t=0.35 

T/t=0.5 

T/1=1.0 



1.0 2.0 3.0 4.0 



FIG. 3. Optical conductivity (a) in the x A -y~ 
dered low-temperature phase for T = 0.05,0.1 and 0.2i and 
(b) in the high-temperature regime with gapless orbital ex- 
citations. Parameters as in Fig. 1 and V = O.lt. The inset 
in (a) shows the Drude weight D c , the sum rule S and the 
incoherent spectral weight Su versus temperature. Inset (b): 
Optical conductivity in the band model for the hole concen- 
tration x = 0.6 for various temperatures (P = 0.02t). 
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Figure 3a shows the regular part <?(u>) at low tem- 
perature in the orbital ordered phase (the Drude peak 
~ D c S(lo) is not shown). The conductivity is given in 
dimensionless form po<r(oj), where po = fiaje 2 A 
typical value is po ~ 0.2 • 10 -3 ficm. The three data sets 
for T — O.Obt, O.lt and 0.2t almost coincide, consistent 
with the saturation of orbital order below T — 0.2t (Fig. 
1). The finite frequency absorption a{uj) for x = 0.6 in 
Fig. 3a shows a well developed gap uj opt ~ 2t, which 
is due to the pronounced gap in the orbital excitation 
spectra. We note that in the x 2 -y 2 ordered phase at low 
doping studied in Ref. jl(J there is also a pseudogap in 
er(<j), which is however smaller in size and only observ- 
able at low temperatures, because of the smaller orbital 
excitation gap. The inset in Fig. 3a shows the tempera- 
ture dependence of the optical sum rule S — — ^ (H*™), 
the absorption at finite frequency S u = -^i a(ui)dui 
and the Drude weight D c . The latter quantity is deter- 
mined from the sum rule S = D c + ]l0| . The weight 
of the Drude peak strongly increases below T = 0.5i and 
saturates below T — 0.2t, i.e. similar to the evolution of 
orbital order in Fig. 1. At low temperature about 60% 
of the spectral weight is in the Drude peak, which may 
however depend somewhat on the cluster size. 

There is a dramatic shift of spectral weight in <j(u>) 
at higher temperatures when entering the orbital liquid 
phase (Fig. 3b). The conductivity is characterized by a 
broad peak centered at u = 0, rcminescent of a conven- 
tional Drude peak with a natural broadening T ~ 0.5i. 
We attribute this intrinsic broadening to the finite den- 
sity of low-energy orbital excitations in the orbital liquid 
phase (Fig. 2b). There is still some structure in <j(u>) 
near lu ~ 2t close to the onset of absorption in the or- 
dered low-T phase. 

The conductivity for the planar model calculated using 
the band picture is shown in Fig. 3b (inset). In this ap- 
proach the finite frequency absorption is due to interband 
scattering, i.e. a consequence of the off-diagonal hopping 
t^. At low temperature the results for <t(uj) obtained by 
neglecting correlation effects are not very different from 
the exact diagonalization result. Both methods yield a 
clear gap and the main absorption is close tow ~ 2t. The 
onset of absorption in the band picture is near u) ~ It, 
i.e. about 0.5£ lower than in the cluster calculation. The 
temperature dependence in the band picture is, however, 
quite distinct. The filling in of spectral weight at low 
energy occurs only gradually, i.e. even at temperatures 
as high as T = It there is only little low-energy spectral 
weight. The temperature scale of the transition from the 
x 2 -y 2 ordered phase to the orbital liquid is a consequence 
of strong correlations, and not contained in the band pic- 
ture, where the natural energy scale is much larger, i.e. 
given by the band width W = 6t. 

Planar x 2 -y 2 correlations are favored by the kinetic en- 
ergy of holes, since the partially filled x 2 -y 2 band has the 



largest overlap. The motion is quasi-2D since the hop- 
ping along the c-direction vanishes by symmetry. More- 
over the antiferromagnetic coupling between layers in the 
A-phase in combination with double exchange would pre- 
vent coherent c-axis motion. 

In a cubic system there are now two possibilities: (i) 
Since x 2 -y 2 correlations are not preferred compared to 
y 2 -z 2 and z 2 -x 2 correlations in the corresponding per- 
pendicular planes, the system remains in a disordered 
quantum liquid state. This is the orbital liquid state 
proposed by Nagaosa and coworkers [0. (ii) Another 
possibility is spontaneous symmetry breaking selecting 
one particular set of planes, e.g. the formation of x 2 -y 2 
long-range order. We note that in the planar model the 
cubic symmetry is explicitely broken, therefore case (i) 
does not occur for this geometry. In the planar model 
the orbital liquid develops at higher temperature due to 
thermal fluctuations. The orbital ordering temperature 
T°° increases with doping. 

We believe that in cubic manganites both cases are re- 
alized in different doping regimes. At moderate doping 
concentration the orbital liquid state appears to be the 
proper ground state in combination with the isotropic 
ferromagnetic spin structure. This is consistent with the 
broad absorption in <j{u>) observed by Okimoto et al. and 
others. These spectra can be compared with our result 
for the orbital liquid in the planar model (Fig. 3b) which 
is stabilized by thermal fluctuations. 

The observation of the A-phase at high doping concen- 
tration by Kawano et al. |l| and Akimoto et al. |!| seems 
to be a realization of case (ii) where symmetry is broken 
spontaneously. Our prediction for this phase are gapped 
orbital excitations with energy ~ 3xt. The optical con- 
ductivity of this quasi-2D metallic state is expected to 
consist of a large Drude peak, while the finite frequency 
absorption shows a pronounced gap due to the gapped 
orbital excitations. 
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